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We study the single site SU{N) lattice gauge theory with Nf = 2 adjoint Wilson fermions for 
values of up to 53. We determine the phase diagram of the theory as a function of the hopping 
parameter K and the inverse 't Hooft coupling b, searching for the region in which the center 
symmetry is unbroken. In this region the theory is equivalent to the infinite volume theory when 
N ^ °°. We find a region of values of K on both sides of for which the symmetry is unbroken, 
including both light physical quarks and masses ^ i?{\/a). This is surrounded by a region with 
a complicated sequence of partially broken phases. We calculate Wilson loop expectation values 
and find that using < 53 it is possible to extract the heavy-quark potential at small distances 
(1-3 links) but not at longer distances. For this, larger values of A^, or lattices with more sites, are 
needed. 
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1. Introduction 

A key property of\aige-N lattice gauge theories is volume reduction, first introduced by Eguchi 
and Kawai (EK) [1]. They showed that, under some assumptions, Yang-Mills theory on an infinite 
lattice satisfies the same Dyson-Schwinger equations as the theory reduced to a single space-time 
point. The crucial assumption was that the center symmetry of the action cannot be spontaneously 
broken. This assumption was soon shown to be invalid [2]. Over the years, many alternatives have 
been proposed to restore EK equivalence, but until recently, none has been fully satisfactory. ' 

The recent revival of this topic was triggered by the work of Refs. [4, 5] which showed that 
volume reduction was an example of broad class of\aige-N orbifold equivalences, holding both in 
the continuum and on the lattice. If one can construct orbifold projections between the parent and 
daughter theories (in our case the large-volume and single-site theories, respectively) then the two 
become equivalent forN^ °°, assuming certain conditions hold. The most non-trivial of these con- 
ditions is, as above, the preservation of the center symmetry. The authors of Refs. [5, 6] proposed 
two possible fixes that might stabilize the ground state to preserve the center symmetry. They calcu- 
lated the perturbative potential for the eigenvalues of the Polyakov loop and showed that it becomes 
repulsive if one adds Nf > \/2 massless Dirac fermions in the adjoint representation (with peri- 
odic boundary conditions).^ In perturbation theory this leads to a center-symmetric ground state, 
with eigenvalues spread uniformly ai^ound the unit circle. Apart from being interesting in their own 
right, theories with Nf adjoint fermions are connected by a chain of large-A'^ equivalences to QCD 
with 2Nf fundamental flavours. This opens up the possibility of using single-site simulations to 
give non-perturbative insight into the strong interactions (up to 1 /N corrections). 

In the following we focus on a single-site SU (N) lattice model with adjoint Wilson fermions- 
the Adjoint Eguchi-Kawai or AEK model. The partition function is 

^ = J D[i7, VA,i|/]f'^^^=""^==+^>='"^^^*'^^\ (1.1) 
where 5gauge is the single-site equivalent of the Wilson action: 

5gauge = 2Nb £ ReTrU^UyUj^Ul (1.2) 

{b = -k^ being the inverse 't Hooft coupling) and Dw is the single-site Wilson Dirac operator: 
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(1.3) 



Using massless fermions is crucial in the perturbative analysis of Ref. [6] but a lattice sim- 
ulation by two of us of the Nf = I AEK model found that, for A'^ < 15, center symmetry is sta- 
bilized not only by light fermions but also by fermions with mass I /a [8]. This is much larger 
than the expected upper limit of \/{aN) [9, 10]. If this continues to hold for N ^ °°, then the 



'For an overview see, e.g., Ref. [3]. One highly successful alternative approach due to Narayanan and Neuberger is 
the reduction to a small volume of size ^ (Ifm)^ rather than to a single site. 

^The other possibility is the trace-deformed EK model, which, however, becomes very complex when one wants to 
compactify more than one space-time dimension (see e.g. Ref. [7]). 
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Nf = \ AEK model provides a realization of the original Eguchi-Kawai idea. This result was 
confirmed, and given a semi-analytic understanding (using arguments going beyond perturbation 
theory) in Ref. [11]. Simulations using massless overlap fermions have also found evidence for 
reduction [12]. There has also been work on the Nf = 2 theory also suggesting that reduction 
holds [13]. 

2. Phase diagram otNf = 2 Adjoint Eguchi-Kawai model 

From now on we focus on the AEK model with Nf = 2. In addition to the above-mentioned 
connection to QCD with 4 fundamental flavors, the large volume theory to which it might be 
equivalent is of phenomenological interest as a theory which might exhibit conformal, or near- 
conformal, behavior in the IR. The N = 2 version of this theory has been extensively studied (see 
Ref. [14] and references therein). 

We use the Hybrid Monte Cai^lo algorithm adapted to work on a single-site lattice (for details 
see Ref. [15]). The CPU time scales approximately as A^^ for heavy fermions and A'^^-^ for light 
fermions. ^ This is a significant improvement over the ^ N^ growth for the Metropolis algorithm 
used in Ref. [8], and allows us both to reach higher values of A'^ and reduce statistical eiTors. 

To analyze the phase diagram we have performed scans in the K — b plane using values of A'^ up 
to 30. In addition, we have done several high-statistics runs at selected points in parameter space 
using A'^ < 53. We have performed simulations with values of Z? up to = 200 but we have mostly 
focused on the region b e [0, 1]. In particular, the two values we have looked at with most care are: 

• Zj = 0.35, which is close to the value used in typical lattice simulations (for N = 3 we have 

j8 = 6/^2 ^ 6 3)^ 

• b = l.O, which is a fairly weak coupling (j8 = 18 for A'^ = 3) at which one should be able to 
compare the results to perturbation theory. 

Apart from analyzing the plaquette, we have used two types of observables to study center 
symmetry breaking: 

1. General "open loops": 

K„ = ^Tr [^'^ [/f U'^' , with «^ = 0, ±1 , ±2, . . . (2. 1) 

where U^" = U^". These are the general order parameters for the breaking of the center 
symmetry. The simplest example of such loops are the Polyakov loops =TrU^. 

2. Eigenvalues of link matrices: the single-site gauge transformation — )• D.U^D.^ obviously 
leaves the set of eigenvalues unchanged. In the center symmetric -phase one expects the 
distribution of phases of the eigenvalues to be invariant under translations by 2nn/N. 

■'The contributions to the scaling are A'"' for fundamental link matrix multiplication (adjoint link matrices are never 
explicitly constructed in the updating procedure), ~ A' for the growth in the number of molecular dynamics steps for 
constant acceptance, and, for light fermions, ~ A'"-^ for the increase in the number of CG iterations. 
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Figure 1: Sketch of the phase diagram for the Nf = 2 AEK model. The positions of phase boundaries are 
approximate, and depend somewhat on A^. The funnel gets narrower very slowly when moving to larger 
values of b. We have also found evidence for a Z3 phase on the left-hand side of the funnel for b>5. 

At weak coupling one also expects that the links are close to being simultaneously diago- 
nalizable. We thus apply the gauge transformation that diagonalizes one link and look at the 
elements of the other links. Indeed, we find that for > 1 the diagonal elements of the links 
dominate in this basis while as we move to stronger couplings the off-diagonal elements 
grow the links become increasingly non-commutative. 

The resulting phase diagram is pictured in Figure 1 . Note that this diagram is qualitatively 
similar to the one from Ref. [8] for the Nf = 1 AEK model. The K" = line corresponds to the 
original Eguchi-Kawai model with broken center symmetry at weak coupling. As the fermions 
become lighter a "funnel" of center-symmetric phase appears on both sides of the critical region. 
This is the region in which volume reduction holds. Around it there are regions in which the 
center symmetry is broken — moving away from the funnel the breaking is to increasingly small 
subgroups until we reach a completely broken phase. We label the phases by the approximate 
remnant of the center symmetry. In particular, a "Z," phase has / clumps in the histograms of the 
phases of link eigenvalues. In the broken phases there are strong correlations between links in all 
directions — this is why we label these phases with Z,- instead of Z^. For a more detailed analysis 
of the phase structure, see Ref. [15]. 

An important issue is how the width of the funnel depends on A'^ and b. As shown in Fig. 1, we 
denote the lower limit of the funnel by Kf, and the corresponding quark mass by iiif = 1/ (2Ky) — 
1/(2k"c)). We first discuss the b dependence. The analysis of Ref. [11] predicts that at weak 
coupling the center-symmetric and completely broken phases are separated by many phases with 
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Figure 2: Scans of the average plaquette for = 23 at = 0.35 and 1. 



partial breakings (which is in agreement with our results) and that the funnel closes as my ~ b~^/^ 
for ft — )• oo. We have performed simulations for = 10 up to ft = 200 and our results are consistent 
with this form. It is important to note, however, that, although the funnel closes, it does so in such a 
way that quarks of arbitrarily high physical mass, m/ = amphys, remain inside the funnel for ft — )• oo. 

Using runs up to A/^ = 30, we find that Kf increases with N, so that the width of the funnel 
shrinks. It is thus very important to go to even higher N in order to perform a reliable extrapolation 
of Kf to N ^ oo, and thus to establish whether the funnel (at fixed ft) has finite width at N = °° 
or not. At this stage, the strongest evidence we have for a finite width is obtained by studying the 
large-N limit of the plaquette within the putative funnel region. We have done this extrapolation at 
several values of for ft = 0.35 and 1.0, using N up to 53. We find that the extrapolated plaquette 
for K away from is almost independent of K, and furthermore hes close to the value obtained in 
the large-A'^ pure gauge theory [15]. This is what we would expect if reduction holds, because the 
heavy quarks lead only to a small renormalization of the plaquette. 

Another interesting aspect of the phase diagram is the order of the critical line. Large volume 
SU{2) simulations with adjoint fermions find a second order transition line emanating from ft = oo, 
jc = 1/8 which becomes first order below ft « 1/4. This structure is related to the claimed presence 
of an IR fixed point. Although N = 2 is far from the N = o° limit, it is expected for the adjoint 
theory that its properties should depend only weakly on A'^. 

In fact, we find a very different phase diagram, in which there is a clear first-order transition 
for ft up to at least 1. This is illustrated by the results shown in Fig. 2. The discontinuity in the 
plaquette drops as A'^ increases, but, for ft < 1, we find that it remains non-zero in the N ^ oo 
limit. For larger values of ft our results are less definitive. Although we still see a discontinuity, its 
magnitude decreases rapidly with increasing ft, and it is difficult to determine whether it remains 
as N ^ °°. Nevertheless, it is clear that the first-order transition extends to much larger ft than in 
the large-volume N = 2 simulations. The simplest interpretation of our results is that the theory is 
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Figure 3: Log-linear plot of 1 x L Wilson loop versus L for L <N. 



confining in the IR, with lattice artifacts leading to a first-order transition at Kc which extends to 
b ^ °°. 



3. Towards physical quantities 

The ultimate goal of using volume-reduced simulations is to calculate physical quantities that 
one can apply to large volume theories. One such quantity is the heavy quark potential, which can 
be obtained using rectangular Wilson loops. In the single-site theory these become 

WiLuLj) = £ (lTr[/^'[/,^^C/;^'[/f^). (3.1) 

For large L2 we expect W{Li,L2) ~ e^^^^'^^^ but at finite we must keep Lj <N to avoid finite 
effects. 

Figure 3 shows a log-linear plot of 1 x L Wilson loops for several values of A'^. We see conver- 
gence to a common envelope for small L followed by a slow rise. The latter is a finite-A'^ effect — as 
N ^ °° the exponentially falling envelope will extend to L = 00. We stress that this finite -A'^ ef- 
fect is statistically significant. From the minima of the curves, one can estimate how L^ax scales 
with A'^. The dependence is approximately logarithmic, so that calculation of the potential poses 
a significant numerical challenge. The problem becomes worse as we move to larger sepai^ations 
and we are only able to see convergence to the linear envelope (and thus extract the potential) for 
distances of 1-3 links. For larger distances one needs to either use lai^ger values of A'^ or move to 
larger lattices (e.g. 2^). 



4. Summary & outlook 

The single-site lattice gauge theory with 2 flavours of adjoint Wilson-Dirac fermions exhibits 
a rich phase structure, one that is qualitatively explained by the semi-analytic model of Ref. [11]. 
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Our most important result is that, for N up to 53, we find a broad funnel of center-symmetric phase 
where laige-N reduction holds. It extends from strong to very weak coupling and contains fermions 
with masses up to ~ 1/a. The N ^ oo extrapolations of the plaquette suggest that the funnel does 
not close in the laige-N limit but we cannot completely rule out this possibility at present. A 
detailed analysis using even higher is desirable. 

We have also analyzed the order of the critical line Kc, finding a first-order phase transition 
for all values b. This is different from the large volume SU (2) simulations and favors a confining 
scenario over the conformal one for the large-A^ theory. 

To use the large-A^ volume reduction as a practical tool one needs to be able calculate long- 
distance observables. We have presented a calculation of the heavy-quark potential and showed 
that l/N effects limit the precision that one can achieve. 

This analysis can be extended by calculating more physical quantities such as particle masses, 
as well as extending the range over which the heavy-quark potential can be extracted. For the latter 
it is desirable to reduce the "1 /N noise" e.g. by using slightly larger lattices or using twisted model 
as in Ref. [1 1]. In addition, making a careful analysis of the vicinity of the critical line, perhaps 
using improved actions, would be of considerable interest. 
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